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Abstract. The paper deals with the statistical analysis of several data sets as- 
sociated with shape invariant models with different translation, height and scaling 
parameters. We propose to estimate these parameters together with the common 
shape function. Our approach extends the recent work of Bercu and Fraysse to 
multivariate shape invariant models. We propose a very efficient Robbins-Monro 
procedure for the estimation of the translation parameters and we use these es- 
timates in order to evaluate scale parameters. The main pattern is estimated by 
a weighted Nadaraya- Watson recursive estimator. We provide almost sure con- 
vergence and asymptotic normality for all estimators. Finally, we illustrate the 
convergence of our estimation procedure on simulated data as well as on real EGG 
data. 



Statistic analysis of models with periodic data is a mathematical field of great 
interest. Indeed, a detailed analysis of such models enables us to have a satisfactory 
approximation of real life phenomena. In particular, SEMOR models |12] are often 
used to describe a large number of phenomena as Meteorology [23] , road traffic [2], 
[7\ or children growth [8j. Here, we choose to focus our attention on a particular 
class of these models called shape invariant models, introduced by Lawton et al. 
|14) . Periodic shape invariant models are semiparametric regression models with 
an unknown periodic shape function. In this paper, we consider several data sets 
associated with a common shape function and differing from each other by three 
parameters, a translation, a height and a scale. Formally, we are interested in the 
following shape invariant model 



where 1 < j < p and 1 < i < n, the common shape function / is periodic and the 
variables are random, independent and of the same law. The classical approaches 
to estimate the different parameters of the model are to minimize the least-squares 
or to maximize the likelihood of the model when the law of Ei^ is known. Here, we 
propose a new recursive estimation procedure which requires very few assumptions 
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and is really easy to implement. The case where the Xi are equi-distributed deter- 
ministic variables has been considered in [7j, [21j or [23j. When p = 1, ai = 1 and 
Wi = 0, Bercu and Fraysse [T] propose a recursive method to estimate the transla- 
tion parameter 9i. In this paper, we significantly extend their results as we are now 
able to estimate, whatever the value of the dimension parameter p is, the height 
parameter f , the translation parameter 9 and the scale parameter a, respectively 
given by 

[1.2) 

Our first goal is to estimate the translation parameter 6. Estimation of shifts has 
lots of similarities with curve registration and alignment problems jl9]. Analysis 
of EGG curves [21] or the study of traffic data j2], [7] fall into this framework. In 
[8], Gasser and Kneip propose to estimate the shifts by aligning the maxima of 
the curves, their position being estimated by the zeros of the derivative of a kernel 
estimate. In the case where = 1 and Vj = 0, Gamboa et al. [7\ provide a 
semiparametric method for the estimation of the shifts. They use a Discrete Fourier 





Transform to transport the model (1.1) into the Fourier domain. The important 
contribution of Vimond [23j generalizes this study, adding the estimation of scale 
and height parameters. When the parameter 6 is supposed to be random, Gastillo 
and Loubes |2] provide sharp estimators of 6, following the approach of Dalalyan 
et al. [4j in the Gaussian white noise framework. Then, they recover the unknown 
density of 6 using a kernel density estimator. In this work, for the estimation of 6, we 
propose to make use of a multidimensional version of the Robbins-Monro algorithm 
|20] . Assume that one can find a function : — )■ M^, free of the parameter 6, such 
that (j){6) = 0. Then, it is possible to estimate 6 by the Robbins-Monro algorithm 

(1-3) 9n+l = 9n + InTn+l 

where (7^) is a positive sequence of real numbers decreasing towards zero and (T„) is 
a sequence of random vectors such that E[T„+i|J-'„] = 0(^ri) where J^n stands for the 
cr-algebra of the events occurring up to time n. Under standard conditions on the 
function (f) and on the sequence (7^), it is well-known [H] that 9^ tends to 9 almost 
surely. The asymptotic normality of 6'„ may be found in [17] whereas the quadratic 
strong law and the law of iterated logarithm are established in [16]. Results for 
randomly truncated version of the Robbins-Monro algorithm are given in pJiJ . 

Our second goal concerns the estimation of the scale parameter a. The estimation 
of the scale parameters, and more particularly of the sign of the scale parameters, 
is very important as we can see in the numerical illustrations on daily average 
temperatures in [23j. Here, we obtain a strongly consistent estimate of the scale 
parameters a using the prior estimation of 9. 

The last main theoretical part of the paper is referred to the nonparametric es- 
timation of the common shape function /. A wide range of literature is available 
on nonparametric estimation of a regression function. We refer the reader to [5], 
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for two excellent books on density and regression function estimation. Here, we 
focus our attention on the Nadaraya- Watson estimator [15j [25j of /. More precisely, 
we propose to make use of a weighted recursive Nadaraya- Watson estimator jS] of / 
which takes into account the previous estimation of f , 9 and a, respectively by Vn, 
On and a„. It is given, for all x G M, by 



with 

/raj (2^) 

and 



1 y:=iW.,a^){y,,-v,^,,] 



Wnjix) = —K' 



where O^-i^j and j are respectively the j-th component of 6'„_i and a„. 
Moreover, the kernel is a chosen probability density function and the bandwidth 
{hn) is a sequence of positive real numbers decreasing to zero. The main difficulty 
arising here is that we have to deal with the additional term 6n inside the kernel K. 

The paper is organized as follows. Section 2 presents the model and the hypothesis 
which are necessary to carry out our statistical analysis. Section 3 is devoted to the 
parametric estimation of the vector v, while Section 4 deals with our Robbins-Monro 
procedure for the parametric estimation of 6. Section 5 concerns the parametric es- 
timation of the vector a. In these three sections, we establish the almost sure conver- 
gence ofvn, On and as well as their asymptotic normality. A quadratic strong law 
is also provided for these three estimates. Section 6 deals with the nonparametric 
estimation of /. Under standard regularity assumptions on the kernel K, we prove 
the almost sure pointwise convergence of /„ to /. In addition, we also establish the 
asymptotic normality of /„. Section 7 contains some numerical experiments on sim- 
ulated and real data, illustrating the performances of our semiparametric estimation 
procedure. The proofs of the parametric results are given is Section 8, while those 
concerning the nonparametric results are postponed to Section 9. Finally, Section 



10 is devoted to the identifiability constraints associated with the model (1.1). 



2. MODEL AND HYPOTHESIS 

We shall focus our attention on the model given, for 1 < j < p, with p > 2, and 
for all 1 < 2 < n, by (|lT} 

Yij = ajf{Xi - 9j) + Vj + Eij, 

where, for all 1 < j < p, aj 7^ 0. For all 1 < i < n, the noise (^jj) is a sequence 
of independent random variables with mean zero and variances E [efj] = ctj, and 
independent of the random points Xj. In addition, as in [Ij, we make the following 
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hypothesis. 

(Hi) The observation times (Xj) are independent and identically distributed 

with probability density function positive on its support [—1/2; 1/2]. In 
addition, g is continuous on M, twice differentiable with bounded derivatives. 

(7^2) The shape function / is symmetric, bounded, periodic with period 1. 

Our goal is to estimate the parameters a, 9, v as well as the common shape 



function /. However, the shape invariant model (1.1) is not always identifiable. 
Indeed, for a given vector of parameters (a, 6',f) and a given shape function /, one 
can find another vector of parameters (a*,6'*,f*) and an other shape function /* 
such that for all 1 < j < p and for all x G M, 

(2.1) a,/(x - 6,) + V, = a*nx - 9*) + v*. 

Hardle and Marron [IQl or Kneip and Engel [Hj were among the first to discuss 
about identifiability. However, in most papers dealing with shape invariant models, 
the identifiability is not so clear. Indeed, the identifiability conditions proposed by 
some authors could be not as restrictive as it should be. For example, under the 
conditions provided by Vimond |23j, it is possible to find two different vectors of 



parameters (a, 9, v) and (a*, 9*,v*) satisfying their conditions and (2.1 ), while under 



the conditions provided by Wang et al. [24], it is possible to find two different shape 



functions / and /* satisfying their conditions and (2.1). That is the reason why we 
have chosen carefully the identifiability constraints. More precisely, we impose the 
following conditions. 

(-Hs) fix)dx = 0, 

J-1/2 

{Hi) ai = I, 9i = Q and max \9A < 1/4. 

Hypothesis ("Hs) allows us to define uniquely the Vj and the second constraint in 
(7^4) define uniquely the 9j, while ai = 1, 6'i = implies that the first curve is taken 
as a reference. These conditions are well-adapted to our framework. Note that ("Hs) 
and {Ha) could be replaced by (^3) and {H'^ defined as 

.1/2 

("^3) / f{x)dx = and sup |/(a;)| = 1, 

J-1/2 xe[0;l] 

(7^4) 9i = 0, min aj > and max 16*^1 < 1/2. 

i<i<p i<j<p 

An other alternative could be to substitute the hypothesis {H4,) by (H'l) or (H'^') 
defined as 

{H'D ai = 1, 6'i = 0, min aj > 0, max \9j\ < 1/2, 
(H'i') ai >0,9i= 0, a]=v and max \9j\ < 1/4. 

Finally, note that the hypothesis of symmetry of / is not necessary to ensure 



identifiability of the model ( 1.1 ) but this hypothesis makes the estimation of 9 easier, 
as we shall see in Section 4 below. All those identifiability conditions are discussed 
in Section 10. 
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3. ESTIMATION OF THE HEIGHT PARAMETERS 
Via {I-L2) and (^^3), it is not difficult to see that 

E 



Y 



Then, a natural consistent estimator Vn of v is given, for all 1 < j < p, by 



(3.1) 



n 

-E 



In order to establish the asymptotic behaviour of Vn, we denote by Y the vector 
(3.2) 



Y = 



where 



Y-j = ajf{X -9j)+Vj+ej, 

with X a random variable sharing the same distribution as the sequence (X„) and 
for 1 < j < p, Sj sharing the same distribution as the sequence (sij). 

The asymptotic results for Vn are as follows. 

Theorem 3.1. Assume that (T-Li) to (T-La) hold. Then, we have the a.s. convergence 

(3.3) lim Vn — V a.s. 

and the asymptotic normality 

(3.4) V^{vn-v) -^^^pio,^{v)), 

where T{v) stands for the covariance matrix given by 



r{v) = Gov 



Y 



MX) J- 

In addition, we also have the quadratic strong law 



(3.5) 



1 

hm - — — (vi — 

n-*-+oolog(n) 



V f,; - V 



r(v) a.s. 



4. ESTIMATION OF THE TRANSLATION PARAMETERS 

In all the sequel, we introduce an auxiliary function defined for all i e M^, by 

/a,f{X - 9i) 
(4.1) 0(i)=E D{X,t) : 

\apfiX - Op) 
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where D{X,t) stands for the diagonal square matrix of order p defined by 
(4.2) D{X, t) = -i-diagf sin(27r(X - ti)), . . . , sin(27r(X - t^))) . 

Using the symmetry of /, the same calculations as in [Ij lead, for all 1 < j < to 
"sin(27r(X - tj)) 



(4.3) E 



9{X) 



-a, fix - e, 



ajfi sin(27r(0j — tj)) 



where /i stands for the first Fourier coefficient of / 

fi= cos{2nx)f{x)dx. 

J-l/2 

Consequently, 

/aisin(27r(^i -ti))^ 

(4.4) 0(t) = /i : 

Vop sin(27r (6'p - tp))^ 

From now and for all the following, in order to avoid tedious discussion, we assume 
that fi 7^ 0. Roughly speaking, we are going to implement our Robbins-Monro 
procedure as in [Ij for each component of 6. More precisely, for all 1 < j < p, 
by noting (f)j{t) = aj/i sin(27r(6'j - tj)), if \tj - 6j\ < 1/2, {tj - Oj)(j)j{t) < if 
sign(aj/i) > and {tj — Oj)(f)j{t) > otherwise. Moreover, denote K = [—1/4; 1/4]. 
Then, we define the projection of a; G M on by 

{X if |x| < 1/4, 
1/4 if X > 1/4, 
-1/4 if X < -1/4. 
Let (7„) be a decreasing sequence of positive real numbers satisfying 

oo oo 

(4.5) ^7„ = +oo and ^7^<+oo. 

n=l n=l 

For the sake of clarity, we shall make use of 7„ = 1/n. Then, for 1 < j < p, we 
estimate 6j via the sequence (Onj) defined, for all n > 1, by 

(4.6) dn+l ,j — T^K (On,j + 7n+isign (flj/l) T„,+i^ 
where the initial value ^ ^^'i the random vector Tn+i is given by 

: 

The almost sure convergence for the estimator 6n is as follows. 
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Theorem 4.1. Assume that to {Ha) hold. Then, On converges almost surely to 
9 as n goes to +oo. In addition, for I < j < p, the number of times that the random 
variable 6n,j + '-fn+isign{ajfi) Tn+ij goes outside of K is almost surely finite. 

Remark 4.1. At first sight, the estimation procedure needs the knowledge of the 
sign of ajfi for all I < j <p. However, it is possible to do without. Indeed, denote 
by (Onj) the sequence defined, for n > 1, as 

where 

: 

and by the sequence defined, for n > 1, as 
where 

: 

Then, two events are possible. More precisely, for I < j < p, 

= ^3 lim 1^- .| = 1/4 a.s. 

or 

lim 9~ , = 9j and lim |^^ | = 1/4 a.s. 

Hence, for n large enough, the vector 6^ which is considered is the vector whose 
absolute value of the j-th component is given by min(|^^^|, Nevertheless, for 

the sake of clarity, we shall do as if the sign of ajfi is known. 

In order to establish the asymptotic normality of it is necessary to introduce 
a second auxiliary function ip defined, for all t G MP, as 

(4.8) cp{t) = E[Vit)Vitf] 

where V{t) is given by 

V{t) = diag(sign(ai/i), . . . , sign(ap/i) )L)(X, 



with Y given by (3.2l). As soon as 47r|/i| min |a,| > 1, denote for all 1 < A;, / < p, 
i<i<p 
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Theorem 4.2. Assume that (Hi) to (^4) hold. In addition, suppose that (cij) has 
a finite moment of order > 2 and that 

47r|/i| mill \aj\ > 1. 
i<i<p 

Then, we have the asymptotic normality 

(4.9) V^i9r,-9)^Afp{0,J:{0))- 

Theorem 4.3. Assume that (Hi) to (^4) hold. In addition, suppose that {cij) has 
a finite moment of order > 2 and that 

47r|/i| mill \aj\ > 1. 
i<i<p 

Then, we have the law of iterated logarithm, given, for all w G MP, by 

1/2 / \ 1/2 



lim sup — — w {6n-0) = - lim inf — — w (6'.„ - 

n^oo v21oglogny n^oo \2iogiognJ 

(4.10) = ^/uFE{e)^ a.s. 
In particular, 

(4.11) limsup (—^^ ] (dn - 9) (On - eY = S(^) a.s. 

n-,00 \2\og\ognJ V /V / 

In addition, we also have the quadratic strong law 

1 " ^ 

(4.12) lim- y"{ei-9){9i-9f = ^{9) a.s. 

n-!>oo log n ^—^^ 

Remark 4.2. In the particular case where 

47r|/i| mill \aj \ = 1, 

i<j<p 

it is also possible to show, thanks to Theorem 2.2.12 page 52 of \6], that 

.[^^AOn-9)^mo,m)- 

V log(^) 

Asymptotic results can also be established when 



< 47r|/i| mill |a,| < 1. 



However, we have chosen to focus our attention on the more attractive case 

47r|/i| mill \aj\ > 1. 

Remark 4.3. We clearly have via differentiation 

D(j){t) = — 27r/idiag^ai cos(27r(^^i — ti)), . . . ,ap cos(27r(^p — tp))^ . 
Consequently, the value 

D(j){9) = -27r/idiag(^ai, . . . ,apj 
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does not depend upon the unknown parameter 6 . On the one hand, if the first Fourier 
coefficient fi of f and the vector a are known, it is possible to provide, via a slight 
modification of (4.6), an asymptotically efficient estimator 6 n of 6. More precisely, 
for all 1 < j < p, it is only necessary to replace 'jn = ^/n in (4.6) hy 7„ = "^j/n 
where 

1 

Then, we deduce from the part 10.2.2 of the hook of Kushner and Yin [13] page 331 
that On is an asymptotically efficient estimator of with 

C 



(4.13) 

where for all 1 < k,l < p. 



n{On - 0) 



M\(),l{0) 
47r2|/i|2|afcai|' 



with (p{0) given by (4.8). On the other hand, if fi and a are unknown, it is also 
possible to provide by the same procedure an asymptotically efficient estimator On of 
replacing fi by its natural estimate 



(4.14) 



fl,r. 



n 

-E 

n ^ 



cos(27rXj) 



and replacing the value by the estimate a„j' which is defined in the next Section. 

Remark 4.4. As in it is also possible to get rid of the symmetry assumption on 
f . However, it requires the knowledge of the first Fourier coefficients of f 



fi 



1/2 



cos(27ra;)/(x) dx 



and 



91 



-1/2 



1/2 



sin(27ra;)/(x) dx. 



-1/2 



On the one hand, it is necessary to assume that fi^O or gi ^ 0, and to replace in 
aatrix D 

A(X,t) 



(4.1) the diagonal matrix D{X,t) defined in (4.2), by 

1 



(4.15) 



9iX) 



diag 5(X,ti 



six,t. 



where, for all I < j < p, 

6{X, tj) = fi sin(27r(X - tj)) - gi cos(27r(X - tj)). 



Then, Theorem 4-1 is true for the projected Robbins-Monro algorithm defined, for 
all 1 < j < P; by 

On+ij = iiK[On,j + sign(aj)7„+iT„+ij ), 



where the initial value Oq G and the random vector T„+i is given by 

^^n+1,1^ 

Tn+1 = A{Xn+l,On] 



.Y 



n+l,p/ 
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On the other hand, we also have to replace the second function ip defined in (4.8) by 

^{t) = E[W{t)W{tf], 

where W{t) is given by 

W{t) = diag(sign(ai/i), . . . , sign(ap/i)) A(X, t)r 
Then, as soon as 47r(/^ + gf) min \aj\ > 1, Theorems 



i<i<P 



4.2 



and 



4.3 hold with S(^) 



kl 



given for all 1 < k,l < p, by 

27r{ff + gl){\a,\ + \ai\)-l 

In the rest of the paper, we shall not go in that direction as our strategy is to make 
very few assumptions on the Fourier coefficients of f . 

5. ESTIMATION OF THE SCALE PARAMETERS 
Henceforth, we introduce an other auxiliary function defined for all t G M^, by 

/ai/(X - 9,) 

(5.1) ipit) = E C{X,t) : 

\aJ{X - 9p) 

where C {X, t) is the diagonal matrix of order p, given by 

(5.2) C(X,t) = ^diag('cos(27r(X-ti)),...,cos(27r(X-tp))). 



As for (4.4), we have 



(5.3) 



m = h 



Then, it is clear from Theorem 



4.1 



^ai cos(27r(6'i — tijY 



\flp coii{2-n{9p — tp))^ 
that il){9n) tends to il){9) = fia. Hence, denote 



by (an) and (an) the sequences defined, for > 1 and for all 1 < j < p, by 



(5.4) 

and 
(5.5) 



1 " 

777 X] 



cos(27r(Xi - 9i^ij)) 



nfi ^ giXi 



n 



cos(27r(Xj - 6*^-1^)) 



Y 



where fi^n is given by (4.14). Let Ip the identity matrix of order Ci the first 
Euclidian vector of W and Mp the square matrix given by 

(5.6) Mp = Ip- aef . 

Then, the asymptotic behaviours of a„ and of follows. 
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Theorem 5.1. Assume that (Hi) to {T-La) hold. Then, we have 

(5.7) lim a„ = a a.s. 

n— >-+oo 

and 

(5.8) lim a„ = a a.s. 

and the asymptotic normalities 

(5.9) y^(a„-a) ^A/;(0,r(a)), 
and 

(5.10) (H„ -a)-^Mp (0, Mpr(a)Mj) , 



(5.11) r(a) = ^Cov {C{X, e)Y) . 



where T[a) stands for the covariance matrix given by 

1 

Jl 

In addition, we also have the quadratic strong laws 

n 

(5.12) lim - — (ai — a)(ai — a) =T(a) a.s. 

n-5>+oolog(?2) ^-—^ 

1 " _ 

(5.13) lim - — — — y (oj — a) (oj — a) = MpT(a)MT a.s. 

n-i>+oolog(n) 



6. ESTIMATION OF THE REGRESSION FUNCTION 

In this section, we are interested in the nonparametric estimation of the regression 
function / via a recursive Nadaraya- Watson estimator. On the one hand, we add 
the following standard hypothesis. 

(T^s) The regression function / is Lipschitz. 



and we suppose that /i is known (see Remark 6.1 below). 

On the other hand, we follow the same approach as in [T]. Moreover, for more 
accuracy, we consider a weighted version of the Nadaraya- Watson estimator 



(6.1) fn{x) = ^UJj{x)fn. 

j = l 

where, for all 1 < j < p, 

(6.2) ijjj{x) = ijjj{—x), ojj{x) > and ^^a;j(x) = 1 

i=i 

1 V" (W- ■('r'\ -^ W- (V: . 

(6.3) fn,j{x) 



<3 YJ!=i{w^,,{x) + w,,{-x)) 

and 
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(6.4) 



hr 



-K 



hr. 



The bandwidth (hn) is a sequence of positive real numbers, decreasing to zero, such 
that nhn tends to infinity. For the sake of simphcity, we propose to make use of 
hn = l/n" with a G ]0, 1[. Moreover, we shall assume in all the sequel that the kernel 
i^T is a nonnegative symmetric function, bounded with compact support, satisfying 

/ K{x) dx = 1 and / K'^(x) dx = u'^ . 

Jr Jr 

Theorem 6.1. Assume that (Hi) to (^5) hold and that the sequence {eij) has a 
finite moment of order > 2. Then, for any x G [—1/2; 1/2], we have 



(6.5) 



lim fn{x) = f{x) 



a.s. 



Theorem 6.2. Assume that (Hi) to (Hs) hold and that the sequence {Sij) has a 
finite moment of order > 2. Then, as soon as the bandwidth (hn) satisfies hn = 1/n" 
with a > 1/3, we have for any x G [—1/2; 1/2] with x ^ 0, the pointwise asymptotic 
normality 



(6.6) V^nUnix) - f(x)) -^MlO, 



E 



a]uj]{x) 



^^i^^'j + g{9j - x)) 



In addition, for a; = 0, 
(6.7) v^(/n(0) - /(O)) 



AT 0, 



a 



,=1 (^j9[^v 

Remark 6.1. If fi is unknown, it is necessary to replace the estimator fn{x) defined 
in (6.1) by fn{x) given by 



fn{x) = ^ujj{x)f, 



n,3 



X 



where fn,j{x) is defined, for all x G [—1/2; 1/2] and for all n > 1, by 



fn,j{x) 



6.1 



Etim,ix) + w,,i-x)) 

and 



6.2 are still true for fn{x) with the same 



One can observe that Theorems 
asymptotic variances. 

Remark 6.2. The choice of the weights u)j{x) could be important. Intuitively, the 
asymptotic variances in (6.6) and (6.7) are minimal if for all 1 < j < p, UJj{x) is 
inversely proportional to the variance of the noise. More precisely, the Lagrange 
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Multiplier Theorem gives us the values of the weights for which the asymptotic vari- 



ances in (6.6) and (6.7) are minimal under the constraint (6.2). They are given, 



for all 1 < j < p and for all x E [— 1/2; 1/2] , by 



Uj{x) 



where 



mj{x) 



T.l=imk{x) 
a] {g{9j + x) + g{9, - x)) 



For these values, the asymptotic variances in (6.6) and (6.7) are respectively given, 
for X ^ 0, by 

-1 



and for x = 0, by 



1 + a 



^m^ix) 



E 



mj(0) 



1 + a 2 



7. SIMULATIONS 

In this section, we illustrate the asymptotic behavior of our estimates on simulated 
data as well as on real data. 



7.1. Simulated data. We consider data simulated according to the model (1.1) 



where 1 < j < p and 1 < 
sen the height parameters v 



< n with p = 5 and n = 2 000. We have cho- 
(0, 1/3, —1, 2, —9/10) , the translation parameters 



6 = (0, 1/5, -1/20, -1/7, 1/6) and the scale parameters a = (1, -4, 3, -5/2, -2)' . 
In addition, the noise (Sij) is a sequence of i.i.d. random variables with A/" (0,1) 
distribution. The random variables (Xj) are simulated according to the uniform 
distribution on [—1/2; 1/2] and the regression function, whose fi = 1/2, is given for 
all X G [-1/2; 1/2] by 

5 



/(x) = cos(2fc7rx) 



fc=i 

The simulated data are given in Figure [ij The results of the estimates of the vectors 
V, 6 and a by 6n and are given in Figure [2} The true values are drawn in 
abscissa and their estimates in ordinate. One can observe that the true values and 
their estimates are very close, showing that our parametric estimation procedure 
performs pretty well on the simulated data. 



Moreover, using convergences (3.4), (4.9) and (5.9), one can obtain confidence 
regions for the parameters v, 9 and a. For instance, for n = 2 000, if we denote by 
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Figure 1. Simulated data 




Figure 2. Estimation of v, 6 and a 

In{v2)) In{(^z) and ln{0"i) the Confidence intervals of f2, 9?^ and 04, for a risk 5%, one 
have precisely 

In{y2) = [-0.1160; 0.4452], 
4(^3) = [-0.1211; 0.0226], 
4(04) = [-2.7225; -2.1224]. 
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The lengths of these intervals are respectively 0.5612, 0.1437 and 0.6001. Conse- 
quently, the lengths of In{v2), Ini^s) and /n(a4) are small, which confirm the good 
performance of our parametric estimation procedure. All these confidence intervals 
are drawn in Figure [3j 




200 400 600 aOO 1000 1200 1400 1600 1800 2000 



Figure 3. Confidence intervals of f2, ^3 and 04 



For the estimation of the regression function /, we have chosen a = 9/10 for the 
bandwidth sequence (/i„). Moreover, the kernel K considered is the uniform kernel 



on [—1; 1] and for all 1 < J < p. 



Uj {X ) 



1/p. The estimation of the regression 
and the estimation of / by fn,i 



function / by /„ is given on the left side of Figure 
is given on the right side. Furthermore, it follows 
that for n = 2 000 and for all x G [—1/2; 1/2], a confidence interval for f{x) is given 
by 



Tom convergences (6.6) and (6.7) 



Kjx) 



fn{x) - QfS ^-^ , fn{x) + Qp- 



whereas one can deduce from convergences (3.3) and (3.4) of Theorem 3.2 of ^LJ that 
for n = 2 000 and for all x G [—1/2; 1/2], a confidence interval for f{x) is given by 



■Jn{x) 



7- , \ Vn{x,9n,l) :r / N , 'Vn{x,9n,l) 
fn,l{x) -qi3 Jn,l{x)+q/3 
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where qp stands for the quantile of order < /3 < 1 of the A/'(0, 1) distribution and 
{x,9n) and v^{x,9n,i) are respectively a consistent estimator of the asymptotic 



variance w'^{x,9) in Theorem 6.2 and of v'^{x,9i) in Theorem 3.2 of [T]. In our 
particular case, z/^ = 1/2, and a numerical calculation leads to 



w\x,9) 



' 0.0114 


if 


- 1/2 < a; < -23/50 and 23/50 < a; < 1/2, 


0.0108 


if 


- 23/50 <x< -9/25 and 9/25 < a; < 23/50, 


0.0099 


if 


-9/25 <x < -17/50 and 17/50 < x < 9/25, 


' 0.0086 


if 


- 17/50 < a; < -31/100 and 31/100 <x< 17/50, 


0.0083 


if 


- 31/100 < a; <0 and < a; < 31/100, 


0.0166 


if X 


= 0, 



and 



v\xM) 




if x^O, 
if x = 0. 



Roughly speaking, for all x G [—1/2; 1/2], the order of the asymptotic variance 
v'^[x^9\) obtained from the estimate jn,\ is ten times greater than the order of the 
asymptotic variance w^(x, ^) obtained from /„. In addition, in this case, the optimal 
variance given in Remark 6.2 is, for all x G [—1/2; 1/2], of order 10"'^ which is anew 



ten times smaller than if^(x, 9\ The confidence intervals K^ix) and Jjy^^ are drawn 
in red in Figure [s One can observe in Figure [4] that the estimate /„ is closer to the 
function / than jn,\- More precisely, the estimation of / by is better than the one 
by /n,i because the lengths of the confidence intervals Kn{x) are smaller than the 
ones of the confidence intervals J„(x) as one can see in Figure [5| which is due to the 
order of the two variances f ^(x, 9i) and w^{x, 9). In our case, the largest confidence 
interval Kn{x) is of length 0.3460 (for x = 0) and is almost three times smaller than 
the smallest confidence interval Jn(x) which is of length 0.9723 (for x = —0.09 and 
X = 0.09). In particular, this justifies the choice of taking a weighed version of the 
Nadaraya- Watson estimator. 
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Figure 5. Confidence intervals for / 



7.2. Modeling of ECG data. This section is devoted to real EGG data. An EGG 
is a recording of the electrical activity of the heart over a period of time, as detected 
by electrodes attached to the outer surface of the skin. A typical EGG consists of 
a P wave, followed by a QRS complex, and a T wave. Here, we consider two sets 
of EGG data, one corresponding to a healthy heart and the other to a heart having 
arrythmia. These data are extracted from the MIT-BIH Database, and they are 
represented respectively in Figure [6] and in Figure [7j 

For each EGG recording, one consider that the heart cycle of interest, that is to 
say the cycle PQRST, is roughly the same at each beat. One also consider that 
every heart cycle are noised and that the white noise is independent of the typical 
shape we want to estimate. After an appropriate segmentation of the EGGs, one 
observe signals of same length such that each of them contains an unique PQRST 
cycle. The segmentation of the EGGs is done by detecting the maximum of each 
QRS complex and centering the segments around this maxima. It is very important 
to have segments of same length in order to ensure periodicity. A well-adapted 
method for the segmentation is the one proposed by Gasser and Kneip [8j. After 
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Figure 6. ECG of a healthy patient 
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Figure 7. ECG of a patient having arrythmia 



the segmentation of the two ECGs, we obtain p = 18 and every segments of length 
77, = 83 for the healthy heart and p = 15 and n = 91 for the ill heart. Then, our goal 
is to estimate the typical shape of each ECG, corresponding to the function / in the 
model (1.1). Firstly, the estimation for the healthy heart is going to show that our 
model is well adapted for the problem of modelling an ECG signal. Indeed, for the 
healthy heart, a good approximation of the heart cycle is to take the average of the 
p different signals, whereas the different parameters v, 9 and a of the model (1.1) 
are trivial. The estimation of the typical shape of the ECG for the healthy heart by 
our procedure is on the right hand-side of Figure [8} whereas the left-side shows one 
on the original ECG signal. The comparison between the two plots shows that our 
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estimation procedure performs pretty well and that our estimation procedure can 
be useful for the modeling of an ECG signal. 




Figure 8. One signal and its reconstruction 

Secondly, we want to estimate the typical shape of the nonhealthy heart, ploted 
in Figure [7j One see that the electric activity is more irregular than for the healthy 
heart, and a simple averaging may lead to a mean cycle that does not correspond 



to the typical shape of the ECG. More precisely, we suppose that the model (1.1) 

Yij = CLjfiXi -9j)+ Vj + Eij, 
fits the data. The segmentation allows us to have a common shape function / 1- 



periodic. Moreover, / is nonsymmetric, but we already saw in Remark that our 
procedure still works for nonsymetric shape function. The parameters a, 9 and v 
correspond to the deformation due to the arrythmia relative to the common shape 
/ we want to estimate. For more accuracy, we have to choose one of the p curves as 
a reference, that is to say where for one 1 < j* < a^. = 1, 9j* = and Vj* = 0. 
For this choice, we consider a criterion of residual variance. More precisely, we first 



consider the model (1.1 ) with Oi = 1, 6'i = and Vi = 0. From this model, we apply 
our procedure to estimate the different parameters a, 9 and v and the shape function 
/ respectively by a„, 9n, Vn and /„. With these estimates, we then calculate the 



vector whose j-th component (a^ i)j is defined by 



1 " 

^ i=l 



Then, we make use of the same procedure by changing the curve of reference, and 
we finally obtain p vectors of length p 



^2 ^2 



Finally, the choice of the curve of reference for the modeling of the ECG signal is 
given by taking 



j* = arg min | 

i<j<P 
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where ||.||i corresponds to the /^-norm. Therefore, we model the ECG by 

Yij = ajf{Xi -Oj) + Vj + Eij 

where 

1 < i < n, 1 < j < p and aj* = 1, 6j* = 0, Vj* = 0. 
On our data set, the implementation of this method shows that j* = 3 and 

= 0.6095. 

The result for the estimation of the typical shape / is given in Figure [9j On the right 
side of Figure [9] one can compare the estimation of / obtained with our estimate 
and with the simple average signal. One can observe that our estimation procedure 
is better because the P wave is well estimated. 




Figure 9. Estimation of / by /„ (in red) and with the mean average 
signal (in black) 



8. PROOFS OF THE PARAMETRIC RESULTS 



8.1. Proof of Theorem 3.1, The convergences (3.3) and (3.4) follow from the 



standard law of large numbers and the standard central limit theorem for martingales 



with independent increments. Moreover, we immediately deduce (3.5) from Theorem 
2.1 of [3]. □ 



8.2. Proof of Theorem 4.1 The result follows from Theorem 2.1 of [1]. 



□ 



8.3. Proof of Theorem 4.2 Our goal is to apply Theorem 2.1 page 330 of Kush- 
ner and Yin [13j. First of all, as 7„ = 1/n, the condition on the decreasing step 
is satisfied. Moreover, we already saw that 6n converges almost surely to 6. Con- 
sequently, all the local assumptions of Theorem 2.1 of Kushner and Yin |13| are 



satisfied. In addition, it is not hard to see from (4.7) that 



a.s. 
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Moreover, the function (f) is continuously differentiable. Hence, (f){9) = and D(j){9) 
is the square diagonal matrix defined by 

D(l){e) = -27r/idiag (ai, . . . , Op) . 

By noting /„ the identity matrix, the condition 47r|/i| min |a, | > 1 implies that 

i<i<p 



Dm + T^ip 



1 

2' 

is a negative-definite matrix. Furthermore, we have that, for all 1 < fc,/ < p, 

E [sign{akfi)Tn+i,ksign{aifi)Tn+i,i\J^n] = <^{On)k,i a.s. 
which leads to 

lim E [sign(afc/i)T„+i,fcsign(ai/i)T„+i,;|j:„] = (6*)^; a.s. 
Consequently, if we are able to prove that the sequence (Wn) given by 

\\6n — 
Wn = — 

In 

is tight, then we shall deduce from Theorem 2.1 of |13) that 



where for all 1 < A;, / < p, 

= mk,i r^exp ((1 - 27r|M(|a,| + |a,|))t) dt = ^ru!^^^^2'\ h T 
Jo ^ ' ^7r|/i|(|afc| + \ai\) - 1 

Therefore, it remains to prove the tightness of the sequence {Wn). Let (V^) be the 
sequence defined, for all n > 1, by 

(8.1) Vn=Wn-e\\\ 

and the sequence defined, for all 1 < J < J5, by 

(8.2) T^ = sign(a,/i)T.,,-. 
Then, we clearly have 

Ki+l = ll^n.+ l — 

= I It^Kp (On + Jn+lT^+l^ - T^Kv {0) \ \^ 

as T^KP = {t^Ki ■ ■ ■ , t^k)^ is a Lipschitz function. It follows that 

Vn+i < v; + 7'+il|r^+ill' + 27„+i<^„-^,r^+i> a.s. 
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By taking expectation in the previous inequality, we obtain that there exists a 
constant M > such that 

(8.3) E[K+i|-^n] <Vn + il+^M + 27„+i <en-e, E[T;+i| j;] > a.s. 
Moreover, ( |4.7 ) together with ( |8.2 ) lead to 

(8.4) E[T;,+i|J-„] = 5p(a)0(^, 
where 

S:p{a) = diag (sign (ai/i) , . . . , sign (op/i)) 



Hence, we deduce from (8.3) and (8.4) that 



(8.5) ElWn+ilJ'n] <— + in+iM + 2 < 6^,-9, Spia)^ fe) > a.s. 

7n+l V / 

Moreover, a Taylor expansion of (p allows us to write 

<9n-9, Sp{a)^ (9n^ > = <9n-9, 27r/i5'p(a)diag (ai, . . . , ap) (^9 - 9„^ > 
(8.6) + f\<9r,-9, 5'p(a)diag {a,, . . . , a^) V (9^ (9 -9^] >, 



where for all t ^ 9, 

V (t) = diag ^^""^^^ ~ ^'^^ ~ ^""^^^ ~ ^'"^ ^^""^^^ ~ ^^^^ ~ ^""^^^ ~ 



'1 — (-i 



0p 



Moreover, the equality 



where 



/iS'p(a)diag (oi, . . . , a^) = L(a) 



L(a) = diag(|/iai|,...,|/iap|). 



together with (8.6) lead to 



<9n-9,Sp{a)(l)i9j > = -2ni9„ -9) L{a) [9^ 



T 



n-9) L{a)V{9n] {9n-9 



Hence, ( |8.5[ ) can be rewritten as 

(8.7) EiWn+ilJ'n] < {l + in)Wn + lnM-A7r(9n-9Y L{a)(9n-9 



T 



Moreover, as 



2{9n-9) L{a)V{9n] [9 



L{a) > min \ajfi\Ip, 



we deduce from (8.7) that 



(8.8) ElWn+ilJ'n] < Wn + 2q^J¥n + M^n-2[9n-9) L{a) [9^ - 9 



RECURSIVE ESTIMATION IN A CLASS OF MODELS OF DEFORMATION 



23 



where 



2g = 1 — 477 1 /i I min \a 

i<j<p ■ 



which means that g < 0. By the continuity of the function V, one can find < e < 
1/2 such that, if ||t - 6'|| < e, 



2|/i| min \a. 

^<j<p 



rip < V(t) < 0. 



Moreover, let An and Bn be the sets An = {\\6n ^ 0\ \ ^ ^} and 

n 

Bn = fl^A: 



k=m 



with 1 < m < n. Then, it follows from (8.9) that 

(8.10) 0<-2|/i|min|a,|V(^„)U. <- 

i<i<p 

where Ip is the identity matrix. Hence, it follows from the conjunction of (8.8) and 

(8.10) that for all n > m, 

(8.11) < Wn\BA^ + qin)+lnM. 

Since -Bn+i = Bn H Bn+i C -B„, and we obtain by taking the expectation on 

both sides of (8.11) that for all n > m, 

(8.12) E[Wn+llB„^,] < (1 + q-fnMWnlBj + InM. 

Finally, following the same lines as in the proof of Theorem 2.2 in [Ij, we obtain 
that for all ^ > 0, it exists K > {] such that for m large enough, 

sup ¥{Wn > K)<i 

ri>m 



which implies the tightness of (Wn) and completes the proof of Theorem 4.2 □ 



8.4. Proof of Theorem 4.3[ The vectorial law of iterated logarithm given by 
(4.10) follows from Theorem 1 of [16], while the vectorial quadratic strong law given 
by (4.12) can be obtained from Theorem 3 of [16j. □ 



8.5. Proof of Theorem 5.1 Recall that for all 1 < j < p. 



and 



Then, it is clear that 



1 ^ cos(27r(X,-g,_i,,)) ^^ 



1 COs(27r(X, - g,_ij)) 

1 " 

in = — Y ^ ("^'' 

'^•'^ i=l 
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1 

^fl,n i=l ^ ^ 



.1^ 



We also have the decompositions 
(8.13) 
and 



ttn - a = -^Sn{a) + -^Rn{a), 



(8.14) an -a 

with 

and the remainder 
Moreover, 



1 



nfi,n 



Sn{a) + (/i - /i,n) a) H —Rn{a 

^ ^ ^ nfi,n 



Sn{a) = Y,iC{X^,e)Y,-^a) 



1=1 



n 

Rn{a) = J2 (C*!^- ^-i) - C (X„ d)^ Y,. 

1=1 



n , 

{fi-fi,n)a = J2[f^ 
1=1 ^ 



cos(27rXj) 



Yli a 



-eJSn{a)a 



Hence, we deduce from (|8.14|) that 
(8.15) 



ttn - a = —l—MpSn{a) H h—Rn{a) 

nfi,n nfi^n 



where the matrix Mp is given by (5.6). In addition, for all 1 < j < 
(8.16) 
where 



Ac,: 



1,31 



and 



Firstly, since 



Acij = cos 27r(Xi - 6'i_ij) - cos {2Tx{Xi - Oj)) . 



E[C(X,,^)r,|J-,;_i] = /ia, 
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the sequence {Sn{a)) is a vectorial martingale with independent increments. For all 
n > 1, its predictable variation (S'(a))„ is given by 

n 

(^(a))„ = J2^[{C{X„e)Yi-f,a){C{X„e)Yi~haf\J^,_,] 

i=l 
n 

= J]cov(c(x„0)r,|j-,, 



i-i 



1=1 



Then, it is clear that 



lim 



(5(a)), 



T(a) a.s. 



where T{a) is given by (5.11). 



Secondly, since 



E 



1/2 



-1/2 



COS ( 27r(a; — Oi_ij) ) dx — cos (27r(x — 6'^)) dx 



1/2 



1/2 



0, 



K i-l 



the sequence (i?^ ^(a)) is a square integrable martingale whose predictable variation 
is given by 

Then, as cos is a Lipschitz function, we have 

|Aci,j| < l^j-ij ~ ^il- 

Consequently, since g does not vanish on [—1/2; 1/2], there exists a constant C > 
such that 

r Ac^ . 1 

(8-17) E < C - 



^7^(^. 



Then, it follows from ( 4.12[ ) together with the previous inequality (8.17) that 
(8.18) (i?2(a))„ = 0(log(n)) a.s. 

Therefore, we deduce from the strong law of large numbers for martingales given 
e.g. by Theorem 1.3.15 of [6j that, for all I < j < p, 



(8.19) 



Moreover, (i?^^(a)) is also a square integrable martingale whose predictable varia- 
tion is given by 

(i?|(a))„ = a|(i?|(a))„. 
Then, we immediately deduce from (8.18) that 

(8.20) (/?3(a))„ = 0(log(n)) a.s., 

and from the strong law of large numbers for martingales that, for all 1 < j < p, 

(8.21) RiAa) = o{\og{n)) a.s. 
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Afterwards, for all 1 < j < p, with the change of variables u = x — 6j, 

cos(27r(x — 9i-ij)) — cos {2n{x — 6j)) ) f{x — 6j)dx 



E 



-1/2 

1/2-e 

-1/2-8 



cos(27r(u + — 9i-ij)) — cos (2™) ] f{u)du 



Then, the elementary trigonometric equality 

cos(27r(M + 6'j — 6'j_ij)) = cos(27ru) cos(27r(6'j — Oi^i j)) — sin(27rM) sin(27r(^j — 6'j_ij)), 
the symmetry and the periodicity of the function / lead to 



E 



Ac,; 



fiX,-e^)\T,_, = /i cos(27r(e, - - 1 a.s 



Moreover, for all |a;| < 1/2, we have 

I cos(27rx) - 1| < 27r^x^ 

which implies that 



(8.22) 



IE 



r Ac 



/(X,-^,)|J-,_i \ <2nU9,-9,^,A a.s. 



Moreover, we have the decomposition 



Rljia) = An.iia) + Bn.iia 



where 



and 



S„,,(a) = 5^E -^/(X,-^,)|J-,_ 



It follows once again from the quadratic strong law (4.12) together with (8.22) that, 
for all 1 < j < p, 

(8.23) Bnj{a) = O (login)) a.s. 

Moreover, for all 1 < j < p, {Anj{a)) is a square integrable martingale whose 
predictable variation {Aj{a))n satisfy 

" r Ac^ 1 

i=i ^yy^i) 

As the shape function / is bounded, we deduce from (4.12) together with (8.17) that 

{A,{a))n = (login)) a.s. 

Therefore, we can conclude from the strong law of large numbers for martingales 
that, for all 1 < j < p, 

(8.24) A^jia) =0 (login)) a.s. 
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Finally, we infer from ( |8.19D , ( |8.2lD together with ( |8.23D and ( |8.24D that, for all 
1 < j < P, 



(8.25) 



Hence, one obtain from (8.13) that 

1 



(8.26) 



Rnjia) = O (log(n)) a.s. 

log(n 



an - a = —rSn{a) + O 
nfi 



n 



a.s. 



and from (8.15) that 
(8.27) 



a„ — a 



Mpbn[a) + (J 



nfi,n 



n 



a.s. 



Consequently, as /i „ converges almost surely to /i, (5.7) and (5.8) follow from the 



law of large numbers for martingales and (5.9) and (5.10) follow from the central 



limit theorem for martingales and Slutsky's lemma, while one can obtain (5.12) 



(5.13) from Theorem 2.1 of [3j 



□ 



9. PROOFS OF THE NONPARAMETRIC RESULTS 



9.1. Proof of Theorem 6.1. For x G [—1/2; 1/2], denote by fn,j{x) the sequence 
defined for n > 1 and 1 < j < p, by 



(9.1) 

We can rewrite 

(9.2) 

where 

and 



fn,j{pC^ an,j fn,j i^-^) ■ 



fn,j{x) = fliix) + flAx), 



Er=i(w^M(^)+w^M(-^)) 



^ Er=i(w^M(^)+w^M(-^)) 

On the one hand, it follows from Theorem 3.1 of [Ij that for any x G [—1/2; 1/2], 

lim flAx) = ajf{x) a.s. 

On the other hand, the almost sure convergence of Vi-ij to vj as i goes to infinity 
implies by Toeplitz lemma, that for any xG[— l/2;l/2], 

lim flAx) = a.s. 

Hence, one can conclude that 

(9.3) lim fn,j{x) = ajf{x) a.s. 
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Consequently, as a^j converges almost surely to aj 7^ as goes to infinity, it 
follows that 

(9.4) 



lim fnjix) = f{x) a.s. 



Finally, (6.1) with (9.4) allow us to conclude the proof of Theorem 6.1 



□ 



9.2. Proof of Theorem 6.2 , We shall now proceed to the proof of the asymptotic 
normality of /„. We have, for all x G [—1/2; 1/2], 



fn{x) -fix) = ^UJj{x) (7nj(x) - f (x) 

(9.5) 
where 



i=i 

p 



Qn,j{x) = ttnj idnjix) +gn,j{-x)) 

Mn,j{x) = Mn,j{x) + Mn,j{-X), 
"^njix) Pn,j{x) ~\- Pn j(y 3^), 

Qn,j{x) Qn,j{x) ~\~ Qn,j(^ ''')) 

"^nji^x) Rnji^x) -\- Rn,j{^ x)^ 

^n,j{x) = Sn,j{x) -\- Snj{ 3;), 



with gn,j{x), Mn,j{x), Pn,j{x), Qn,j{x), Rn,j{x) and Sn,j{x) givcu by 

1 " 

9n,jix) = -^Wij{x), 
i=l 

n 



i=l 



n 

Pn,{x) = a„,,5^H^,,(x)(/(X,-^,_i,,)-/(x; 



i=l 

n 



QnM = a„,,^l^„(x) (/(X,-^^,)-/(X,-^,„i 

i=l 

n 

Rn,j{x) = {ttj - an,j) Wij{x)f{Xi - 9. 

i=l 

n 

Sn,j{x) = ^ Wij{x) {Vj - Vi^ij) . 



<3/ 1 



1=1 
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Firstly, (6.28) of [IJ together with the almost sure convergence of a„ to a as n goes 
to infinity, lead to 

(9.6) lim Gn,j{x) = ttj {g{Oj + x) + g{Oj — x)) a.s. 

In addition, we obtain from (6.32) and (6.35) of [Ij that, for a > 1/3, 

(9.7) Vl^{x) = o(ni+") a.s., 

(9.8) Ql/x) = o(ni+") a.s. 

Hence, for all x G [—1/2, 1/2], we find that 

(9.9) lim a/^V u;,(x) ^"'^'^^^ + QnA^) ^ q 
Secondly, as the shape function / is bounded, it follows that 

Rn,j{x) = O ^aj - an,j\ ^ Wij{x)^ a.s. 

Hence, 

n 

= 0{n) a.s., 

i=l 

ensures that 

(9.10) Rnj{x) = O {n\aj — cinjl) a.s. 



In addition, we can deduce from (8.26) that, for all 1 < j < p 



|aj -a„,j| = C I — I a.s. 



which via (9.10), leads to 

(9.11) Rn,j{x) = O (^-^n \og{n)j a.s. 
Consequently, 

(9.12) 'Kjix) = o(n^+") a.s. 
Thirdly, we have the following inequality 

(9.13) |5„,,(x)| < A„,,(x) + S„,,(x) 
where 

n 
i=l 

and 



Sn,i(x) =J2'^i,MWiAx)\J^i^,i 
1=1 
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with Ci j = \vj — We deduce from (6.34) of [Ij together with the Cauchy- 

Schwarz inequahty and the quadratic strong law given by (3.5) that 

(9.14) S„,,(x) = o(v^ fE4.) ) = ^ (\^M^) a.s. 



Moreover, the sequence {Anj{x)) is a martingale whose predictable variation is given 

by 

(A,(x))„ = O |'fj4^.E[Wt,(x)|J-.-i]j a.s. 



Consequently, we obtain one again from the quadratic strong law (3.5) that 



n 



(9.15) (A(a:))„j = O j = ^ ^^sH) 

which allows us to show, from the strong law of large numbers for martingales that 
for any 7 > 0, 

(9.16) Al/x) =o{n" login a.s. 
Therefore, 

Sl.{x) = o (n" log(n)2+T) + O {n \og{n)) 

(9.17) = o(n^+") a.s. 

We are now in position to study the asymptotic behavior of the dominating term 
A4n,j{x). For all x G [—1/2; 1/2] and for all 1 < j < p, the sequence {Ain,j{x)) is a 
square-integrable martingale whose predictable variation is given by 

n 

(A^,(x))„ = [{W,^j{x) + W,^,{-x)f I . 

1=1 

We deduce from (6.37) of [1] that we have, for x 7^ 0, 

(9.18) hm ^-^^^^ = ^ {g{e, + x) + g{e, - x)) . 
and from (6.38) of [1] that, for x = 0, 



.2, ,2 



(9.19) Ito <^iMli.4i21<,(9,). 

Moreover, according to (6.39) of [Ij, as (^jj) has a moment of order > 2, Lindeberg 
condition is satisfied for A^„ j(x). We can conclude from the central limit theorem 
for martingales given e.g. by Corollary 2.1.10 of [6j that for all x G [—1/2; 1/2] with 

(9.20) %M ^ Ar(0, ^ + X) + g{B, - x))) , 
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while, for x = 0, 
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(9.21) 



Finally, it follows from (9.20) and (9.21) and the independence of to- 
gether with the previous convergence (9.6) and Slutsky's theorem that, for all 
X e [-1/2, 1/2] with X ^ 0, 



(9.22) -^^y 

while, for x = 
(9.23) 



Uj{X} 



2 P 

V 1 + 



^^1^1^] +x)+ 9{Gj - a;)) / 



,=1 



vergences (9.22) and (9.23 



Then, the conjunction of (|9.9), (9.12), (9.17) together with the two previous con- 

and Slutsky's theorem let us to conclude that, for all 

a]u][x) 



,2 P 



^^7^x^] (si^j + ^)+ 9{0j - x)) 



X e [-1/2, 1/2] with x^O 

V^(/„(x) - fix)) Ar(o, -2 
while, for x = 0, 

y/nhniw) - fm Ar(o, ^ E ^ 

V 1 + a ^-^ a - 



' ^^M(0) 



which completes the proof of Theorem 6.2 



□ 



10. APPENDIX 

In order to prove every identifiability conditions, let us consider that for a given 
vector of parameters (a, 9, v) satisfying (7^4) and a given shape function / satisfying 
(7/2) and ("Hs), one can find another vector of parameters (a*, 9*, v*) satisfying (^4) 
and an other shape function /* satisfying {'H2) and [T-Lz) such that for all 1 < j < p 
and for all x G M, (2.1) is true, that is to say 

(10.1) a,f{x - 9,) + V, = a*r{x - 9*) + v* 



First of all, the periodicity of / and ("Hs) lead to Vj = v*. Then, (10.1) becomes 

(10.2) a,f{x-9,) = a*r{x-9*). 

For j = 1, the identifiability constraints ai = and 9i = 9^ enables us to show 
that f{x) = f*{x). Then, (10.2) can be rewritten as 

(10.3) a, fix -9,) = a* fix -9*). 
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Hence, by denoting I2 the integral of the square of / 

f'^{x)dx, 







and by squaring and integrating into (10.3), we obtain that 

(10.4) ap,= {a*Yh, 
leading to 

(10.5) a|=(a*)'. 

If aj = a*, then f{x — 6j) = f{x — 9*) and it follows from the constraint 

max \9j\ < 1/4 
i<i<p 

that 6j = 6*. Otherwise, if aj = —a*, then f{x — 6j) + f{x — 6*) = 0, which clearly 
leads to the identitity f{x) = f{x + 2{6j — Oj))- Therefore, the constraint 

max \6j\ < 1/4 
i<i<p 

implies that 6j = 9*. Then, f{x) = 0, which is impossible. Finally, we have shown 
that 

V = V* , 6 = 6*, a = a* and f = f* , 



leading to the identifiability of the model (1.1 ). This reasoning fits for the other sets 
of identifiability constraints. □ 

Acknowledgements. The author thanks Bernard Bercu for all his advices and for 
his thorough readings of the paper. 
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